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Abstract 

We present an analytic study of the finite size effects in Sine-Gordon model, based 
on the semiclassical quantization of an appropriate kink background defined on a 
cylindrical geometry. The quasi-periodic kink is realized as an elliptic function with 
its real period related to the size of the system. The stability equation for the small 
quantum fluctuations around this classical background is of Lame type and the 
corresponding energy eigenvalues are selected inside the allowed bands by imposing 
periodic boundary conditions. We derive analytical expressions for the ground state 
and excited states scaling functions, which provide an explicit description of the flow 
between the IR and UV regimes of the model. Finally, the semiclassical form factors 
and two-point functions of the basic field and of the energy operator are obtained, 
completing the semiclassical quantization of the Sine-Gordon model on the cylinder. 
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1 Introduction 



Quantum field theory on a finite volume is a subject of both theoretical and practical 
interest. It almost invariably enters the extrapolation procedure of numerical simulations, 
limited in general to rather small samples, but it is also intimately related to quantum field 
theory at finite temperature. It is therefore important to increase our ability in treating 
finite size effects by developing efficient analytic means. In the last years, a considerable 
progress has been registered in particular on the study of finite size behaviour of two 
dimensional systems. Also for these models, however, an exact treatment of their finite 
size effects has been obtained only in particular situations, namely when the systems 
are at criticality or if they correspond to integrable field theories. At criticality, in fact, 
methods of finite size scaling and Conformal Field Theory P1I2] permit to determine many 
universal amplitudes and to extract as well useful information on the entire spectrum of 
the transfer matrix. Away from criticality, exact results can be obtained only for those 
integrable theories described by a factorized and elastic scattering matrix [HI Ej which, 
on a finite volume, can be further analysed by means of Thermodynamical Bethe Ansatz 
[5-9]. This technique provides integral equations for the energy levels, mostly solved 
numerically. In all other cases, the control of finite size effects in two dimensional QFT 
has been reached up to now either by conformal perturbation theory or numerical methods 
as, for instance, the one proposed in [11 . 

The aim of this paper is to study the finite size effects of a two dimensional massive 
QFT by using a different approach, i.e. the non-perturbative semiclassical expansion 
formulated in the infinite volume case by Dashen, Hasslacher, Neveu jl2j and by Goldstone 
and Jackiw . Apart from some issues which make such an analysis an interesting 
subject in itself, the main theoretical motivation of this work consists in the possibility of 
obtaining analytic results for the form factors and the energy levels at a finite geometry. 
In integrable cases, this adds to the above techniques (see also [101), whereas for non- 
integrable models it is an efficient alternative to perturbative or numerical studies. As 
a matter of fact, in the infinite volume case, semiclassical methods have proved to be, 
together with Form Factor Perturbation Theory [T^], ideal tools in the analysis of non- 
integrable quantum field theories (see, for instance, Ref. [T7j). 

Form factors at a finite volume of local operators in both integrable and non-integrable 
theories have been studied in one of our previous papers |18j . These quantities enter the 
spectral density representation of correlation functions which need, however, another set 
of data for their complete determination, precisely the energies of the intermediate states 
at a finite volume. This paper is mainly devoted to fill this gap, that is, to face the 
problem of a semiclassical computation of the energies Ei(R) of vacua and excited states 
as functions of the circumference R of a cylindrical geometry. Notice that, isolating a 
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factor 1/R in front of the E^Rj's (simply due to their dimensionality), the remaining 
quantities are scaling functions of the variable r = mR, where m is the lowest mass of 
the considered QFT. 

It is worth to underline an important feature that has come out from the study of the 
semiclassical form factors in infinite volume. As we will discuss later, their accuracy seems 
to extend, somehow, beyond the regime in which they were supposed to be valid. Together 
with the known fast convergency properties of the spectral series and the information that 
can be extracted on energy levels, the above result suggests that the semiclassical method 
may provide a rather precise estimate of finite volume correlation functions, an outcome 
which may be useful for many applications. 

For methodological reasons, we have decided to present the semiclassical computation 
of finite volume energies for a system that admits one of the simplest analysis, the Sine- 
Gordon (SG) model. As we will see, this model is particularly appealing for its simplified 
semiclassical results whereby the significant physical effects we are looking for will not be 
masked by other additional complications. Moreover, due to the integrable nature of this 
theory, its finite size effects have been previously studied by means of Thermodynamical 
Bethe Ansatz E] , and it would be interesting to perform a quantitative comparison 
between these results and the semiclassical ones, in order to directly control their range 
of validity. However, as already pointed out, semiclassical methods apply not only to 
integrable theories and this opens the way to describe analytically the finite size effects 
also in non-integrable models [TH] . 

The paper is organized as follows: in Section El we briefly recall the main ideas and re- 
sults of the semiclassical approach. We also discuss the simplest scaling function in a finite 
volume in order to clarify the nature of divergencies encountered in such computations. 
Section El is devoted to the complete semiclassical analysis of the energies of the quantum 
states in the kink sector of the SG model on a cylinder. In general, this analysis passes 
through the solution of a Schrodinger type equation for a particle in a periodic potential 
and, for the SG model, this corresponds to a Lame differential equation. In this section 
we also discuss how to select the proper eigenvalues inside the band structure of the spec- 
trum in order to determine the energy levels E^R). In Section 0] we compute the form 
factors of local operators by using the semiclassical methods and we comment on their 
properties. Our conclusions and further directions are discussed in Section 03 There are 
also several appendices: Appendix A presents the quantization of a free bosonic theory in 
a finite volume and a comparison of finite-volume and finite-temperature computations of 
the simplest one-point correlation function. Appendix B collects relevant mathematical 
properties of the elliptic functions used in the text whereas Appendix C displays the main 
properties of the Lame equation. 
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2 Semiclassical quantization 



In this section, after recalling the basic equations of the semiclassical quantization, we 
will present the simplest example of a scaling function on a cylindrical geometry, i.e. the 
ground state energy SQ ac (R) of a free massive bosonic field. In a semiclassical quantization, 
SQ ac (R) is the lowest energy level in the vacuum sector of the theory. This example will 
show, in particular, how to handle the divergencies usually encountered in the calculation 
of the scaling functions. 



The main feature of a large class of 2-D field theories with non-linear interaction and 
discrete degenerate minima is that they admit non-perturbative finite-energy classical 
solutions (called kinks or solitons) carrying topological charges Q^ p = ±1. In this paper 
we will concentrate our attention, in particular, on a specific model of this kind, i.e. the 
Sine-Gordon (SG), defined by the potential 



In such theories, the kinks generally interpolate between two next neighbouring minima 
of the potential (vacua) which are constant solutions of the equation of motion (in our 
example 4>sg — ^jf, s — 0, 1), and they exhibit certain particle properties. For instance, 
they are localized and topological stable objects, i.e. they do not decay into mesons with 
Qtop — 0. Moreover, in integrable theories as SG model, their scattering is dispersionless 
and, in the collision processes, they preserve their form simply passing through each other. 

The kinks are static solutions of the equation of motion, i.e. they are time indepen- 
dent in their rest frame, and they can be simply obtained by integrating the first order 
differential equation 



further imposing that 4) c i{x) reaches two different minima of the potential V(<f)) at x — > 
±oo. These boundary conditions, which describe the infinite volume case, require the 
vanishing of the integration constant A. As we will see in the next Section, the kink 
solutions in a finite volume correspond instead to a non-zero value of A, related to the 
size R of the system. 

All the above properties of the kink solutions are an indication that they can sur- 
vive the quantization, giving rise to the quantum states in one-particle sector of the 
corresponding QFT. A direct correspondence among the kink states and the correspond- 
ing classical solution has been established by Goldstone and Jackiw who have shown in 



2.1 DHN method 




(2.1) 
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Ref. ^2] that the matrix element of the field between kink states is given, at leading 
order in the semiclassical limit, by the Fourier transform of the kink background. We will 
discuss this result and its applications in Sect.lU 

At the moment we are mainly concerned with the semiclassical quantization of the 
small fluctuations around kink backgrounds. As it is well known, one cannot apply directly 
to them the standard perturbative methods of quantization around the free field theory 
since the kinks are entirely non-perturbative solutions of the interacting theory. Their 
classical mass, for instance, is usually inversely proportional to the coupling constant (in 
the SG model one has M c i = ^). In an infinite volume, an effective method for the 
semiclassical quantization of such kink solutions (as well as of the vacua ones) has been 
developed in a series of papers by Dashen, Hasslacher and Neveu (DHN) jT2j by using an 
appropriate generalization of the WKB approximation in quantum mechanics (see also 
Refs. ^3] or ^3] for a review). The DHN method consists in initially splitting the field 
<fr(x, t) in terms of its classical static solution of the equation of motion (which can be 
either one of the vacua or the kink configuration) and its quantum fluctuations, i.e. 

<f>(x,t) = <f> d (x) + ri(x,t) , r)(x,t) = ^e^rkix) , (2.3) 

k 

and in further expanding the Hamiltonian of the theory in powers of 77, by keeping only 
the quadratic terms. As a result of this procedure, rj k (x) satisfies the Schrodinger equation 

r] k {x) = uj 2 k r] k {x) , (2.4) 

together with certain boundary conditions. The semiclassical energy levels in each sector 
are then built in terms of the energy of the corresponding classical solution and the 
eigenvalues u)i of the stability equation (|2.4j) . i.e. 

E {ni} = E d + hJ2 Uk + \ ) + 0(h 2 ) , (2.5) 

where n k are non-negative integers. In particular the ground state energy in each sector 
is obtained by choosing all n k = and it is therefore given by 1 

E = E cl + ^Y,^ + ^ ■ ( 2 - 6 ) 

k 

In summary, to each static finite energy solution of ()2.2j) corresponds a tower of quan- 
tum energy eigenstates ()2.5j) representing the O-particle (vacua) and 1-particle (kink), 

^rom now on we will fix ft, = 1. since it is well known that the semiclassical expansion in ft is equivalent 
to the expansion in the interaction coupling. 



d 2 



dx 2 
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and their excitations. The construction of the complete Hilbert space, including the n- 
particle sectors (for n > 2), requires to consider time-dependent multi-kink and breather 
solutions with finite energy. Their semiclassical quantization can be performed with an 
appropriate modification of the DHN method ^2] • 

As we have mentioned in the Introduction, the analytic form of the semiclassical 
scaling functions for the 2-D QFT's admitting static kink solutions can be achieved by 
DHN method suitably adapted to finite size geometry. In the following we will discuss 
in details the results of Sine-Gordon model on a cylindrical geometry which, as we shall 
see, admits the simplest technical analysis. On this geometry — described by a space 
variable compactified on a circle of circumference R and by a time variable t running on 
an infinite interval — the SG model admits quasi-periodic boundary conditions (b.c.) 

2?77r 

<f)(x + R,t) = + — , (2.7) 

where the arbitrary winding number n G Z originates from the invariance of the potential 
(12. lj) under — > + ^ ZE . In particular, since we are interested in the one-kink sector, 
which is defined by n — 1, we will impose the b.c. 

<f){x + R,t) = <f)(x,t) + . (2.8) 

P 

The first step for applying the semiclassical method to this problem is to find the finite 
size analog of the kink solution, satisfying now the b.c.'s (|2.8|) . However, the success in 
constructing the scaling functions depends on whether one is able to solve the correspond- 
ing Schrodinger equation ()2.4|) and to derive an analytical expression for its frequencies u^. 
It turns out that the semiclassical finite size effects in SG model are intrinsically related 
to the simplest (N = 1) Lame equation, which admits a complete analytical study. 



2.2 SG in infinite volume 

The semiclassical quantization of the Sine-Gordon soliton in infinite volume has been 
performed in ^2] • We report here the basic results in order to show how the semiclassical 
technique works in the simplest example and also to introduce the quantities that should 
be obtained in the IR limit of the forthcoming finite volume results. 
The classical (anti)soliton 

<P cl {x) = | arctan ( e ±"*(*-*o)) ( 2 .9) 

is solution of eq. (J2.2J) with A = 0. It connects the two degenerate vacua = and = ^ 
and its classical mass is given by M c \ = 8|p-. Plugging the above expression in ()2.4j) . this 
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equation can be cast in the hypergeometric form by using the variable z = |(l + tanhmx) 
and its solution is expressed as 



rj(x) = izW 1 -^!-^-^ 1 -^ f( 2, -1,1 + ^1-^,^) • ( 2 - 10 ) 



The corresponding spectrum is given by the discrete value uj 2 = (i.e. the zero mode 
associated to the translation invariance of the theory) and by the continuous part uj 2 = 
m 2 (l + q 2 ), characterized by the absence of reflection and by the phase shift 5(q) = 
2 arctan(-). 

The semiclassical correction to the mass is given by the difference between the ground 
state energy of the soliton and the one of the vacuum, with the addition of a mass 
counterterm due to normal ordering of the interaction term in the Hamiltonian: 

2 °° 

M-M cl = l^[ m y/T+ti- + j dx [1 - cos/30 d (x)] , (2.11) 

n „ 



with 



oo 

,2/32 



2 = m 2 p 2 r dk 

8vr 7 ^k 2 + m 2 ' 1 ' ' 



-oo 



The discrete values q n and k n are obtained by defining the system in a large finite volume 
of size R with periodic boundary conditions: 

2nn = k n R = mq n R + 5(q n ) . (2-13) 

Sending R — > oo and computing the integrals, one finally obtains the semiclassical quan- 
tum correction to the mass of the kink 

M = (2.14) 

As it is well known, the exact solution of the quantum Sine-Gordon model shows that 
the coupling constant /3 2 renormalises as 



" 2 ^ = d^- < 2 - i5 > 

Moreover, the exact quantum mass of the soliton is given by M = —, which coincides 
with the above expression (|2.14jl . The equality of the semiclassical and the exact result 
for the soliton mass is a remarkable property of the SG model, on which we will come 
back in the following Sections. 
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2.3 Ground state energy regularization in finite volume 

As shown by eq. (|2.fij) . quantum corrections to energy levels are given by the series on the 
frequencies u n . However, this series is generally divergent (this is the usual UV divergence 
in field theory) and a criterion is needed to regularize it. It is quite instructive to consider 
the simplest example where such divergence occur, i.e. in the calculation of the ground 
state energy ^Q ac (i?) of the vacuum sector of the theory on a cylindrical geometry of 
circumference R. This can be constructed by implementing the DHN procedure for one 
of the constant solutions, for instance 0™ c = 0, imposing periodic boundary conditions 
for the corresponding fluctuations r) vac (x). Obviously, what comes out is nothing else but 
the Casimir energy of a free bosonic field <j)(x, t) with mass m. In this case the frequency 
eigenvalues are fixed to be 

UJ n = V / Pn + ^ 2 , (2-16) 

with p n = 2irn/R and n = 0, ±1, ±2, . . .. 

The ground state energy has to be regularized by subtracting its infinite-volume con- 
tinuous limit: this ensures in fact the proper normalization of this quantity, expressed by 

lim £^°(R) = . (2.17) 

R— >oo 

The ground state energy at a finite volume is therefore defined by 



oo 



Isolating the zero mode, it can be conveniently rewritten as 



(2.18) 



. oo 

oo 



W0 .- + !5^ + £) , -!/*^ + (£)". „ 19) 

n=l Q 

where r = mR. Since the divergence of the series is due to the large n behaviour of the 
first two terms in the expansion 



" 2 + (£> 2 -" + ^)> G?)' < 2 - 20 > 

we begin our calculation by subtracting and adding these divergent terms to it: 



V27T. 

n=l n=l 



n=l 




n=l 
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The first series in the right hand side of the above expression is now convergent, whereas 
the last two terms should be coupled to the analogous ones coming from the integral, 
whose divergencies have to be handled in strict correspondence with those coming from 
the series. Hence, by subtracting and adding the leading divergence to the integral 




(2.22) 



we can combine the last term in this expression with the one in (|2.21|) and implement the 
well known regularization 

= ~ ■ < 2 -») 

However, the first term in ()2.22j) still contains a subleading logarithmic divergence, as it 
can be seen by explicitly computing the integral by using a cut-off A, in the limit A — > oo 



E 

n=0 



n dn 



lim 

a->0 



E 

n=0 



ne 



ne 



dn 



This divergence can be cured by subtracting and adding the term | (^) 2 lnA. By com- 
bining this last term with its analogous in the series we have 



lim [ — 

A->oo \ ^— ' n 

\n=l 



In A = 7s , (2-25) 



where 7# is the Euler-Mascheroni constant, while the remaining part of (|2.24j) with the 
above subtraction is now finite. 

Collecting the above results, the finite expression of the ground state energy on a 
cylinder is then given by 



it r r 
6 + 2 + 4^r 



- (\n — + 7 S - + V (j(2nn) 2 + r 2 - 2irn - -^—] 
Ti V 4tt ' 2 J ^ V 47m/ 

x ' n=l x ' _ 

(2.26) 

It is now easy to see that (J2.26j) fulfills the requirement of modular invariance of the 
theory, which imposes its equality with the TBA expression [Zj 

Sr(R) = , (2.27) 



S 



where 

c(r) = j / ^cosh^ln(l-e- rcoshe ) . (2.28) 

In fact, this integral formula can be expressed in terms of Bessel functions, which admit 
a series representation that directly leads to ()2.26|) (see Ref. jZj). For this theory we obvi- 
ously have c(0) = 1. Moreover, one can also check that the above regularization scheme 
ensures the agreement between the R and L channel calculations of the finite expression 
of the one-point functions (<p 2k ) [21]- The interested reader can find the simplest example 
of these calculations in Appendix 1X1 

Finally, it is worth to note that the result ()2.26|) can also be obtained by using a simpler 
prescription which automatically includes the subtraction of the various divergencies, 
fastening the calculation. This consists in ignoring the divergent part of the integral, 
keeping only its finite part, and in regularizing the divergent series as 



oo 

£ 

n=l 
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oo 

^ n 

n=l 



res 



1 
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(2.29) 
(2.30) 



Formula (J2.29)) is the standard regularization of the Riemann zeta function £(— 1), where 
C( s ) = S^Li ~ > an d usually corresponds to normal ordering with respect to the infinite 
volume vacuum (see, for instance, (221; chapter 4). On the contrary, the regularization 
of the second series is a-priori ambiguous due to its logarithmic divergence, and its finite 
value (|2.3Uj) was chosen according to the above discussion. 



3 Scaling functions on the cylinder 

We will now develop a complete semiclassical scheme to analyse the energy of the quantum 
state in SG model containing one soliton on the cylinder. This can be achieved by applying 
the DHN method to an appropriate kink background. 



3.1 Properties of the periodic kink solution 

In order to identify a kink on the cylinder, we have to look for a static finite energy 
solution of the SG model satisfying the quasi-periodic boundary condition (|2.8|) . For the 
first order equation 

1 / d6ri \ 2 m? , 



+ A) (3.1) 
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a solution which has this property can be found for A > 0. It can be expressed as 

, . , 7r 2 / mix — x Q ) , 9 \ , 9 2 . 

Mx) = ^ + 1 , ° ; ,fc 2 , fc 2 = — — , (3.2) 



P V fc / 
provided the circumference R of the cylinder is identified with 

R = —kK(k 2 ) , (3.3) 

m 

where K(/c 2 ) denotes the complete elliptic integral of the first kind 2 . The parameter xq 
in ()3.2j) represents the kink's center of mass position, and its arbitrariness is due to the 
translational invariance of the theory around the cylinder axis. The behaviour of (|3.2j) as 
a function of the real variable x is shown in Figure ^ 

Ma 



~K(k 2 



_ m.i: 
k 



Figure 1: Solution of eq. (|3.1|) with A > and x = 0. 



The function (|3.2j) has been first proposed in and interpreted as a crystal of 
solitons in the sine-Gordon theory in infinite volume. In our finite volume case, instead, 
(13.21) has to be seen as a single soliton defined on a cylinder of circumference R (given 
by eq. (|3.3j) ). while its quasi-periodic oscillations represent winding around the cylinder. 
As shown in eq. ()3.3|) . there is an explicit relation between the size of the system and 
the integration constant A. It is easy to see that the infinite volume solution (|2.9J) is 
consistently recovered from ()3.2j) in the limit A — > 0, i.e. when R goes to infinity. 

2 The definition and basic properties of K(fc 2 ) and the Jacobi elliptic amplitude am(u, k 2 ) can be found 
in Appendix IbI 
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Sd(R) 



The classical energy of the kink on the cylinder is given by 

E(k 2 ) 



R/2 




J dx 


"l fd(j> d \ 


2 V dx ) 


R/2 





m 



1 — COS , 



k 



+ 



1 - 



k 2 



K(k 2 ) 



(3.4) 

where E(/c 2 ) is the complete elliptic integral of the second kind. In the R — > oo limit 
(which corresponds to k! — > 0, with {k') 2 = 1 — k 2 ), £ c i(R) approaches exponentially the 
correct value = This can be seen expanding E and K for small k' (see Appendix 
B), and expressing the result in terms of mR, which can be itself expanded in k' in virtue 
of the relation (|3.3)l : 



-mR 



Hence the large R expansion of the classical energy is 

32 



£ d (R) = Moo + — 



me 



-mR 



O (e~ 2mR ) . 



(3.5) 



We will comment more on the interpretation of this result in Section | 

Similarly, one can derive the behaviour of £ c i{R) for small r = mR, which corresponds 
to the limit A — >• oo (or k 2 — > 0): 



Sd(R) 



2/T 7T r 

HJ 2 + m J 2 



m 



r v 
27rJ 



7T 



2/3^ 



+ 



(3.6) 



This formula will be relevant in the discussion of the UV properties of the scaling functions 
presented in Sect 13. 31 

Before moving to the quantization of the kink-background (J3.2)) . it is worth to mention 
that another simple kind of elliptic function, which solves eq. (j3.1|) for —2 < A < 0, was 
also proposed in [22] and interpreted as a crystal of solitons and antisolitons in the infinite 
volume SG. This background corresponds as well to a kink on the cylinder geometry but 
satisfying the antiperiodic boundary conditions 



[x + R, t) 



2tt 
J 



<j>{x,t) 



The associated form factors were obtained in ^8]. Although the quantization of this 
second kink solution is technically similar to the one of ()3.2|) presented in the next section, 
it displays however some different interpretative features that justify its discussion in a 
separate publication [TH] . 
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3.2 Semiclassical quantization in finite volume 

The application of the DHN method to the periodic kink ()3.2j) requires the solution of 
eq. ()2.4|) for the quantum fluctuations tj u , which in this case takes the form 

^ + k 2 (iu 2 + l) -2k 2 sn 2 (x,k 2 )^ V Ux) = 0, (3.7) 

where sn(x, k 2 ) is the Jacobi elliptic function defined in Appendix El and we have intro- 
duced the rescaled variables 

mx to 
x = — — , uj = — . (3.8) 
k m 

Due to the periodic properties of (p c i(x) expressed by eq. (|3.2[) . the boundary condition 
(12.8)1 translates in the requirement for t)q(x) 

Vu(% + —r-} = Vu>(%) ■ ( 3 - 9 ) 



Eq. ()3.7j) can be cast in the so-called Lame form, which admits the two linearly indepen- 
dent solutions 

= g( ' + f^ a) , (3.10) 

a[x + iK ) 

where the auxiliary parameter a is defined as a root of the equation 

2 - k 2 

V{a) = — k 2 u 2 . (3.11) 

The Weierstrass functions V(u), ((u) and u{u) are defined in Appendix O where the 
Lame equation and its relation with ()3.7j) are discussed in detail. 

As it is usually the case for a Schrodinger-like equation with periodic potential, the 
spectrum of eq. (J3.7)) has a band structure, determined by the properties of the Floquet 
exponent 

F{a) = 2i[KC(a)-flC(K)] , (3.12) 
which is defined as the phase acquired by r)± a in circling once the cylinder 



V±a (x + 2K) = e ±iF ^r] ±c 
We have two allowed bands for real F(a), i.e. 



x 



0<u 2 <^--l and uj 2 > — , (3.13) 
k l k l 



and two forbidden bands for F(a) complex, i.e. 



u 2 <0 and — - 1 < uj 2 < — . (3.14) 
k z k 
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k 3 



/n UJ 2 



F" 1 ^ F 



F = 3n 
F = 2tt 



TT 



TT 



-- F = 



allowed band 



forbidden band 



allowed band 



forbidden band 



Figure 2: Spectrum of eq. ()3.7j) 



The band < uo 2 < 1 



p^ is described by a = K + iy, where y varies between and K' and, 
correspondingly, F(a) goes from to w. The other allowed band uj 2 > p corresponds 
instead to a = iy and, by varying y, F(a) goes from n to infinity, as it is shown in Fig. |21 
By imposing the periodic boundary conditions (|3.9J) on the fluctuation r)(x), one selects 
the values of uj 2 for which the Floquet exponent is an even multiple of ir, thus making 
the spectrum of eq. ()3.7|) discrete. These eigenvalues are Uq = 0, which is the zero mode 
associated with translational invariance and has multiplicity one, and the infinite series 
of points 



1 



2 - Ar 



V{iy r , 



(3.15) 



with multiplicity two, placed in the band u > p, with y n determined by the equation 



F{iy n ) = 2Ki((iy n ) + 2y n ((K) = 2nn 



n 



1,2, 



(3.16) 



In the IR limit (A — > 0) the above spectrum goes to the one related to the standard 
background ([2.9)1 : the allowed band < uj 2 < p — 1, in fact, shrinks to the eigenvalue 

merges in the continuous part of the 



q — 0, while the other allowed band u 2 > -p 



spectrum Ur — 1 + q . 

It is useful to note that, although the R dependence of the frequencies (13.15)1 is quite 
implicit, since it passes through the inversion of eq. ()3.3jl . nevertheless these are analytic 
functions of R and it is extremely simple to plot them. The corresponding curves, shown 
in Figure El provide an important piece of information, since they are nothing else but 
the energies of the excited states with respect to their ground state Sq(R). 

To complete the analysis, it remains then to determine the finite volume ground state 
energy £q(R) of the kink sector. In analogy with the infinite volume case (see eq. (12.11)1 ). 



13 




this is defined by 

R/2 

oo r. 2 r- 

£o(R) = S d (R) +J2"n(R) J dx[l- cosPM - £™ C (R) . (3.17) 

n=1 -R/2 

Before commenting in detail each of these terms, let's focus first on the main problem 
in deriving a closed expression for Sq(R), which consists in the evaluation of the infinite 
sum on the frequencies uj n (R) or, better, in isolating its finite part. We need therefore 
a method for solving the transcendental equation (J3.16)) for y n (k 2 ) in order to make the 
expression (|3.15|) for the frequencies u n (k 2 ) explicit. As we have already seen for the 
classical energy, two kinds of expansion are possible, one in the elliptic modulus k and 
the other in the complementary modulus k', which are efficient approximation schemes in 
the small and large r regimes, respectively. Here for simplicity we only present the small 
r expansion. By taking into account the series expansion in k for K, Q{u) and V(u) (see 
Appendices iBl and . we are led to look for a solution of eq. (J3.16J) in the form 

oo 

yn (k 2 ) = Y,( k2 yy ( n s) ■ ( 3 - 18 ) 

s=0 

(s) 

Here we give the result for the first few coefficients yn , s — 0, 1, 2: 

Vn ) = arctanh-J- , 
in 

= \y^ ] , (3.19) 

? /2) = 9 ..(0) U 
tin at tin 



64 yn 16(4n 2 - l) 2 
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which are those relevant in the analysis of the UV properties of the scaling function in 
Sect. l3.3l As a consequence, we obtain the following simple expression for the frequencies: 



m 



2n 
~k 



k 2 k* 20n 2 
~ T ~ 64 An 2 - 



9 



0(k e 



Comparing it order by order with the small-/c expansion of eq. ()3.3|) 

r = mR = it k 



we finally obtain the explicit i?-dependence 

U n (R) 27T 

= — n + 



r 

2^ 



71 



+ 



(3.20) 



(3.21) 



(3.22) 



m r \2ir/ An 2 — 1 

It is worth noting that this series expansion in r, which can be easily extended up to 
desired accuracy, efficiently approximates the exact energy levels also for rather large 
values of the scaling variable. Fig.H] shows a numerical comparison between the first 
energy level given by (j3.15j) and its approximate expression ()3.22|) . and for the higher 
levels it is possible to see that the agreement is even better. 

LOx/m 



12 3 



5 6 7 



Figure 4: Comparison between the exact energy level uj\/m given by (J3.15|) (continuous 
line) and the approximate expression (|3.22|) (dotted line). 

With the above analysis, the sum over frequencies in the ground state energy ()3.17|) 
takes the form 

OO / „n _ OO CO 

= — Yn+ — Y — — - + ... 3.23 
m r ^ \2tx) ^ An 2 - I K ' 

n=l n=l n=l 

As we will see below, the subtraction of counterterm and vacuum energy in ()3.17|) leads 
to the cancellation of all the divergencies, producing a finite expression for the ground 
state energy in the kink sector. 
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Moving now to the analysis of the remaining terms in (|3.17|) . a similar series expansion 
can be easily performed on each of them. The classical energy Sd(R), given in eq. (|3.4j) . 
has already be treated in this way in eq. (|H.6|) . The finite volume counterterm (C.T.), 
where the one-loop mass renormalisation is given by 



m 2 (3 2 2tt 



E 



ill R ' / o . (2n7r) 2 

n=-oo \ m z + - 



(3.24) 



and 4> c i is given by (|3.2|) . takes the explicit form 

K (k 2 ) — E (k 2 ) 



C.T. 



m 



k 



E 



v/(2mr) 2 + r s 



The first terms of its expansion in R are then 



C.T. 



1 
1 



r 

4tt 



oo 

E 1 

71=1 



32tt 2 



r 



4 V2tt 



Ef 1 ^ 



n=l 



+ ... 



(3.25) 



(3.26) 



Finally, the vacuum energy £™ C (R) is the one precisely computed in Sec. 12.31 Since its 
role is to cancel certain divergencies present in the other terms of Sq(R), in complete 
analogy with the infinite volume case (see eq. (|2.11jl ). we will now consider its "naive" 
formulation, given by 



£™(R) 



2m^V\^/ +m "2 + r ^ U + 4vr 8 \2tt/ ^ n 3 + " ' 

n=— oo > ' n=l n=l n=l 

(3.27) 

Hence, in the final expression for the ground state energy all the divergent series present 
in the sum over frequencies, in the counterterm and in the vacuum energy cancel out, and 
one obtains 



£ (R) 2ir 7T 1 1 1 
m Vp 2 ' 4 ^2 r_ 8^ 



f \ 3 

2tt. 



^C(3)-i(2 lo g 2-l] 



7T 



2/3 2 



(3.28) 



where we have used jSU] 

oo 

E 



2n 2 



| C(3) - |(2 log2-i; 



8n 3 (4n 2 - 1) 

n=l v ' 

in order to evaluate explicitly the coefficient of the r 3 term. 

Repeating the above calculations, one can also easily write the finite expressions of 
the excited energy levels (|2.5j) . whose series expansion in r is given by 



m 



2tt 
r 



71 

J 2 

2tt. 



n 



1 

4 



/J 2 



|C(3) - ^(2 log 2-1] 



7T 

2/J 2 



8 V2tt 

E fc 



(3.29) 



4n 2 



+ 



16 



where {k n } is a set of integers defining a particular excited state of the kink. 



3.3 UV— IR correspondence 

The semiclassical quantization of the periodic kink (|3.2j) provides us with analytic expres- 
sions, albeit implicit, for the scaling functions in the kink sector for arbitrary values of the 
scale r = mR. These quantities control analytically the interpolation between the Hilbert 
spaces of the ultraviolet (UV) and infrared (IR) limiting theories. It is worth to note that, 
although we obtain them in the framework of a particle-like description proper of the IR 
limit, the kink background (|3.2|) is intrinsically formulated on a finite size, leading to the 
possibility of extracting UV data. Hence, it is important to check whether our scaling 
functions reproduce both the expected results for the IR (r — > oo) and UV (r — > 0) limits. 

Concerning the IR behaviour, we have seen in the previous Sections that in the R — ► oo 
limit all the quantities in exam, i.e. the classical solution, its classical energy and the 
stability frequencies, correctly reach their asymptotic values. In addition, it is also possible 
to perform a simple and interesting analysis of the first correction to the kink mass for 
large R. According to Liischer's analysis [21] , the mass of a particle in a large but finite 
volume has to approach exponentially its asymptotic value, in a way controlled by the 
scattering data of the infinite volume theory. Restricting for simplicity our analysis to the 
leading term in (3 in the kink mass (which, in our approch, is simply given by the classical 
energy), we have then to compare the expansion presented in (|3.5|) with the term that 
dominates Liischer's formula for small j3. This is given by [2U EH] 

M(R) -M 00 = -m b smu k kb R kbk e~ m ^ u ^ R + • ■ • , (3.30) 

where Rkbk is the residue (multiplied by —i) of the kink-breather S'-matrix on the pole at 
9 = iu k kb , i.e. 

Rkbk = -iResS kb (9 = iu\ b ) . (3.31) 

Using the kink-breather S"-matrix |3] 

, sinh 6 + i cos ^ B 2 , 

S kb (e) = . . f , 7 = = ^ 2/a 3.32 

sinned — 2cosy^ 1 — p z /87r 

and selecting its s-channel pole 9* = iu k kh — % (| + ^J, we find 

Rkbk = -2cotg-^ . (3.33) 
lo 

Substituting in f)3.30|) . for small f3 2 we have 

32 

M(R) - Moo = m — e~ r + ■ ■ ■ , (3.34) 







17 



which therefore reproduces eq. (|3.5|) . It is a remarkable fact that the classical energy alone, 
being the leading term in the mass for (3 2 — > 0, contains the IR scattering information 
which controls the large-distance behaviour of £q(R). 

The UV behaviour for r — > of the ground state energy Eq(R) of a given off-critical 
theory is known to be related instead to the Conformal Field Theory (CFT) data (A, A, c) 
of the corresponding critical theory and to the bulk energy term as |5j 

E (R) ~ % (A + A--^-) + BR+--- (3.35) 



R V 12 

where c is the central charge, A + A is the lowest anomalous dimension in a given sector 
of the theory and B the bulk coefficient. 

For the Sine-Gordon model the bulk energy term is given by [HI I2E] 

B = 16— tan J- , 3.36 
7^ 16 

while its UV limit is described by the CFT given by the gaussian action with c = 1 

^cft = \g J d'xd^d^, (3.37) 

where the free bosonic field is compactified on a circle of radius 1Z. The various sectors 
of this CFT are labelled by two integers, s and n: s is the momentum index, while n is 
the winding number, related to the boundary condition imposed on 

<f>(x + R,t) = <f>(x, t) + 2im1l. (3.38) 

Let | s, n) be the states carrying the lowest anomalous dimension in each sector. They 
are created by the vertex operators [21] 

V a , n (z,z) =: exp [ial n <p{z) + ia~ n ^(z)] : , 



i.e. 



where 



s,n) = K,n(0,0) | vac) , (3.39) 



a t,n = ^±^gnK ; 
(x,t) = ip(z)+(p(z) . 



Their conformal dimensions are given by 
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The vacuum sector is described by s = n = 0, with A mc + A vac = 0. If we now define 1Z = 
-^=-p and fix the normalization constant to the value 3 g — 1, then the kink sector in SG, 
defined by the boundary condition ()2.8|) . naturally corresponds to the sector characterized 
by s — 0, n — 1, in which the lowest anomalous dimension is 

Ao,i + A ,i = j 2 . (3.41) 

The conformal vertex operator Vo,i has been put in exact correspondence with the soliton- 
creating operator of SG in Ref. [25] . 

The question to be addressed now is whether the small r expansion of SQ ac (R) and 
£o(R) given by eqs. (j2.2ti|) and ()3.28|) reproduces, in semiclassical approximation, the above 
data controlling the UV limit of SG model. 

For the vacuum sector, comparing (|2.26|) with (|3.35|) . we correctly obtain c = 1 and 
A vac = A vac = 0. We do not expect, however, to obtain the bulk term B relative to 
SG model by looking at ()2.26j) . simply because the semiclassical expression of the ground 
state energy in the vacuum sector applies equally well to any theory which has a quadratic 
expansion near the vacuum state. Namely, apart from the value of the mass m, eq. ()2.26|) 
is a universal expression that does not refer then to SG model. 

The kink scaling function ()3.28|) has instead a richer structure. The obtained scaling 
dimension 

A + A = j 2 (3.42) 

is the expected one for the soliton-creating operator in Sine-Gordon while the central 
charge contribution c = 1 is absent, simply because in ()3.28|) we have subtracted the 
vacuum ground state energy from the kink one 4 . Moreover, the bulk coefficient B = ^ 
present in ()3.28j) correctly reproduces the semiclassical limit of the exact one, given in 
eq. ()3.36|) . In principle, this bulk term should be present in all the energy levels, included 
the ground state energy in the vacuum sector, but its non-perturbative nature makes 
impossible to see it in the semiclassical expansion around the vacuum solution, which is 
in fact purely perturbative. Hence it is not surprising that to extract the bulk energy 
term we have to look at the kink ground state energy, in virtue of the non-perturbative 
nature of the corresponding classical solution. 

Finally, the expression ()3.29j) for the excited energy levels explicitly show their cor- 
respondence with the conformal descendants of the kink ground state. In fact, their 
anomalous dimension is given by 

A {kn) + A {kn} = j 2 +Y,k n n. (3.43) 

n 

3 Note that the usual normalization adopted in the CFT literature is instead g = j-. 
4 The value c = 1 , coming out from the regularization of the leading term of the series on the frequencies 
(13.22(1 . is in fact exactly cancelled by the same term in the vacuum energy. 
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The successful check with known UV and IR asymptotic behaviours confirms the 
ability of the semiclassical results to describe analytically the scaling functions of SG 
model in the one-kink sector. It would be interesting to further test them at arbitrary 
values of r through a numerical comparison with the results of [HJ |§] in an appropriate 
range of parameters. This was not pursued here because the results presently available 
in the literature were obtained for values of (3 which are beyond the semiclassical regime 
and moreover the energy levels were plotted as functions of a different scaling variable, 
i.e. the one defined in terms of the kink mass. We hope however to come back to this 
problem in the future. 

4 Form factors and correlation functions 

The semiclassical scaling functions, derived in Sect. El provide an important information 
about the finite size effects in SG model. As in the infinite volume case, however, the 
complete description of the finite volume QFT requires to find, in addition to the energy 
eigenvalues (|3.29jh the kink form factors and the correlation functions of local operators. 
This section is devoted to the analysis of this problem, i.e. to the determination of the 
finite volume form factors and the corresponding spectral functions. 

4.1 Infinite volume form factors 

It is useful to initially recall some basic definitions and results concerning semiclassical 
form factors for the SG model in infinite volume. As mentioned in Sect. [21 the relation 
between kink solutions and form factors was established by Goldstone and Jackiw [T3] . 
who showed that the matrix element of the field between two asymptotic one-kink 
states is given, at leading order in the semiclassical regime, by the Fourier transform of 
the classical solution describing the kink itself (see also |T3] for further developments). 
This remarkable result, however, had the drawback of being formulated non-covariantly 
in terms of the kink space-momenta. It was refined in fH] with a covariant formulation in 
terms of the rapidity variable 9 of the kink, defined in terms of its energy and momentum 
as E = M cosh 9 , p = Msinh^. In the semiclassical regime, there are moreover further 
simplifications: in fact, the mass of the kink can be approximated by its classical energy 
M ~ M c i = whereas its rapidity can be assumed to be very small, i.e. 9 = 0(/3 2 ), 
thus obtaining E ~ M c i , p ~ M d 9 . Hence, the refined form of Goldstone and Jackiw 
result is given by 

poo 

(9, | 0(0) | 9 2 ) = M d / dae- m ^-^ a (j> d (a) , (4.1) 
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where | Oi) are asymptotic one-kink states. Moreover, it is also possible to prove that the 
form factor of an operator expressible as a function of is given by the Fourier transform 
of the same function of <p c i. For instance, the form factor of the energy density operator e 
can be computed performing the Fourier transform of e c i(x) = ^ 

i^t) +y\u- with this 

covariant formulation, the matrix element (|4.1j) can be continued to the crossed channel, 
and from its pole structure one can easily extract the spectrum of the bound states of the 
theory, even in non-integrable CcLSGS, clS discussed in [13 EE] • 

In ^Hj we have checked that the form factor (j4.1)l obtained from the infinite volume 
kink f)2.9j) reproduces the semiclassical limit of the exact one, derived in [2H]. Here we 
would like to present a more quantitative comparison which permits to conclude that 
formula (|4.1jl . though proven under the semiclassical assumption of small coupling and 
small rapidities, remarkably extends its validity to finite values of the coupling and to 
a large range of the rapidities. Consider, for instance, the form factor of the energy 
operator (up to a normalization N) F{6) = N{ 9 2 | e(0) | 9\ ), whose semiclassical and 
exact expressions are given, respectively, by 

*Uud.(0) = ^ f r ( 4 - 2 ) 



sinh 



4-7T n 



7 cxact(#) = sinh^ G(6) (4.3) 



where £ = 7/871 and 

G(9) = exp 



2 sinh 2 ^ 



50 dt sinh |(l-0 sin 2 || 
t sinh | £ cosh \ sinh t 



(4.4) 



Fig. El shows how, for small values of the coupling, the agreement between the two func- 
tions is very precise for the whole range of the rapidity. Furthermore, the discrepancy 
between exact and semiclassical formulas at larger values of /3 can be simply cured, in our 
example, by substituting the bare coupling (3 2 with its renormalized expression 7 into the 
semiclassical result (|4.2j) . as shown in Fig.H3 Hence we can conclude that the monodromy 
factor (j4.4|) . which is the relevant quantity missing in our approximation, does not play 
a significant role in the quantitative evaluation of the form factor even for certain finite 
values of the coupling 5 . 

As we have already mentioned, the exactness (or very high accuracy, as in this case) 
of the semiclassical results is a peculiar feature of SG model in infinite volume, obtained 

5 It is easy to understand the reason of this conclusion in the above example: at small values of 9 we 
have G{6) ~ 1, whereas for 9 — > 00, when G(9) may contribute, the whole form factor goes anyway to 
zero. Similar conclusion can be reached for all other form factors which vanish at 9 — > 00. 
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F = 0.1 F /3 = 0.5 
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Figure 5: Comparison between the exact function F given by (J4.3|) (continuous line) and 
its semiclassical approximation (|4.2|) (dotted line), at (3 = 0.1 and (3 = 0.5. 



F (3=1 ( a ) F f3=l (b) 
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Figure 6: Comparison, at f3 = 1, between (a) the exact function F given by ()4.3)1 (contin- 
uous line) and its semiclassical approximation (|4.2jl (dotted line), (b) the exact function 
F given by (|4.3|) (continuous line) and its semiclassical approximation ()4.2|) with the 
substitution [3 2 — > 7 (dotted line). 

with the "dressing" j3 2 — > 7. An interesting problem to be studied is whether similar 
phenomena take place for the semiclassical scaling functions and form factors in finite 
volume as well. An indication on this issue could be found by extending to finite volume 
the analysis of higher loop quantum corrections which, in the semiclassical approach, are 
obtained by keeping cubic (and higher) powers of rj in the expansion of V(<f) c \ + 77) |14j . 

4.2 Semiclassical spectral functions on the cylinder 

The generalization of the above construction to the case of finite volume has been proposed 
in [THj, where we have shown how to estimate the leading semiclassical behaviour of the 
spectral function on the cylinder under the same hypotheses of the infinite volume case. 
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In fact, the matrix element 



f(0n) = (Pn 2 \m\P ni ) 



(4.5) 



of the basic field <fi between two kink eigenstates of the finite volume hamiltonian can 
be expressed, at leading order, as the Fourier transform of the corresponding classical 
solution: 



R/2 

f(9 n ) = M(R) J dae lM ^ e " a <p d (a) 

-R/2 



RM(R) 



E 



-iM{R)e n a 



f{0n) ■ 



(4.6) 



(4.7) 



Here we have denoted the states by the so-called "quasi-momentum" variable p n , which 
corresponds to the eigenvalues of the translation operator on the cylinder (even multiples 
of tt/R), and we have defined 8 n as the " quasi-rapidity" of the kink states 



2nn 



Pn 



M{R) sinh# n ~ M{R)9 n 



where M(R) is the classical mass of the finite-volume kink. As shown in |18j . the crossed 
channel form factor can be obtained at leading order via the change of variable 6 — *■ iit — 9: 



F 2 (e n ) = ( o 1 0(0) \p n2 p ni ) = f{m - e n ) , 

and the leading terms in the spectral function on the cylinder are given by 



p{E k , Pk ) = 27i6(E k )5( Pk )\ <0| 0(O)| 0> 



7T 



*(§-2) 



M 2 



Foim 



Pk_ 

M 



(4.9) 



(4.10) 



The procedure described above has been introduced in [18] for the construction of 
form factors for the kink backgrounds in the SG model and the broken 4 field theory, 
both defined on a cylindrical geometry with antiperiodic boundary conditions. In what 
follows we will apply it instead to the case of SG model with periodic boundary conditions. 
The corresponding finite volume form factor ()4.6|) can be written in terms of the soliton 
background (J3.2|) : 



f(0 n ) 



R/2 

M J dae 

-R/2 



iM9„a 



2 

(3 • 



7i 2 ( mx 2 
- am (^ _ ^ _ ; k 



(4.11) 



2?r I M 

J 



R5 



1-5 



M8 n ,0 



Aie n ,o 



cos (M9 n R/2) 



sin (M6 n R/2) 
MQ n R/2 



+ i 



6 n cosh(kK'%6 n ) 
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In order to obtain this result one has to compare the inverse Fourier transform (|4.7|) with 
the expansion |3Ti] 



am m 



71 U 

2K 



1 



i n cosh 



n7r 



sin 



nir- 



u 
K 



(4.12) 



The form factor (j4.11j) has the correct IR limit 6 , and leads to the following expressions 
for F 2 (8) and for the spectral function 7 : 



Am 

Wn 



1 



cosh 



1 " ^„,o 



fc K' — 0, 



cos (M0„ i2/2 



sin (M0 n R/2 
M6 n R/2 



(4.13) 



(1 - <Wo) 



(3 2 (Pn) 



cosh 



fcK' 



+ 



cos (p n R/2) 



sin (p n R/2) 

PnR/2 



(4.14) 



where 6 = in — 6. Note that the finite volume dependence of both the form factor 
()4.13|) and the spectral function (|4.14|) is not restricted to the second term only. The 
k~K'(k 2 ) M(R) factor in the first term carries the main independence, although it is not 
manifest but implicitly defined by eq. (j3.3j) . 

Another quantity of interest is the two-point function ( | e(x)e(0) \ ) of the energy 
density operator. One can calculate it by evaluating its spectral function 

R/2 

p £ {p 2 ) = J dx ( | e(x)e(0) | ) e~ ip - x (4.15) 

-R/2 

in terms of the form factors of e(x), similarly to what we have done for the SG field 
above. In order to find the semiclassical form factor 



}e{0 n ) = (Pnak(O) \p ni ) 

we need to compute the Fourier transform of 



2m 2 



Am 



2 fmx\ 



(4.16) 



(4.17) 



6 The functions cos ( x ^/ 2 ) anc j ^^f^^ can De shown to tend to zero in the distributional sense for 
R — > oo. 

Here we are considering the matrix elements on the antisymmetric combinations of kink and antikink. 



21 



This can be easily obtained from the following expansion 

1 L ff 2 Ancos^ 



sn 2 u 



71=1 K 



and we finally have 



The corresponding semiclassical spectral function is thus given by 

4tt 3 „ ( E n S\ vl 



MEn , Pn) ,_ s ^. 2) _^_. (4 , 0) 

It is worth to mention that it is also possible to obtain the two-point functions of 
certain vertex operators V h (x, t) = e ±t/3b ^ x '^ (for b = |, 1, |, 2, ...), since the required 
Fourier expansion formulas of the type f)4.18|) are known in these cases |3*T] . 



5 Further directions 

In this paper we have shown how the semiclassical methods can provide an analytic 
description of finite size effects in two-dimensional quantum field theories displaying de- 
generate vacua. In particular, we have applied these techniques to study the SG model, 
quantizing its kink solution on the cylinder. The scaling functions of the ground (and 
excited) states, as well as the form factors and two-point functions of different operators, 
allowed us to build the one-kink sector (i.e. Qtop = ±1) of the corresponding Hilbert 
space of states. 

The next step in this program is the extension of the DHN method to describe the 
multi-kink states (Qtop = ±2, ±3,...) as well as the non-vacua ("breather" -like) part 
of the Q t0 p = sector. These states are related to certain time-dependent solutions 
of SG model on the cylinder, i.e. to the finite volume analog of soliton-soliton, soliton- 
antisoliton and breather solutions. Although more complicated from the technical point of 
view, the determination of these classical solutions and the study of their scaling functions 
and form factors is a well stated open problem in the semiclassical framework, which 
deserves further attention. 

One of the advantages of the semiclassical method is that it works equally well for 
both integrable and non-integrable models, if they admit kink-type solutions. In fact, 
we have chosen to test the efficiency of the semiclassical quantization on the example of 
SG model, mainly because it leads to the simplest N = 1 Lame equation. Static elliptic 
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solutions for other models can be easily obtained by integrating equation (2.2) with A ^ 
and appropriate boundary conditions. This was done, for instance in [T^j, where we have 
derived the form factors between kink states in the broken (ft 4 model on the cylinder with 
twisted boundary conditions. In this case, the quantization of the finite volume kink 
involves a Lame equation with N = 2. Lame equations with N > 2 are also expected to 
enter the quantization of other theories. 

Finally, the semiclassical method seems to be suited also for the description of finite 
geometries with boundaries, say with Dirichlet or Neumann boundary conditions. Due to 
the physical significance of this kind of systems, an interesting problem is the semiclassical 
computation of the relative energy levels, a subject that will be discussed in a forthcoming 
publication [THj. 
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A Free theory quantization on a finite geometry 

Let us consider a free bosonic field <fr(x, t) of mass m defined on a cylinder of circumference 
R, i.e. satisfying the periodic boundary condition 

(j)(x + R,t) = (j)(x,t) . (A.l) 

Imposing the equation of motion and the commutation relation 

[<j>{x,t),Il(y,t)} = iS P (x-y) , (A.2) 

where H(x,t) = ^(x,t) is the conjugate momentum of the field whereas 

S P (x) = - ^ X > S P (x + R) = 5 P (x) 

n=— oo 

is the periodic version of the Dirac delta function, we obtain the mode expansion of the 
field <f)(x,t). This is given by 

OO 

(j)(x,t) = E [A n e i{PnX - Unt) + A ] e~ i{PnX - Unt) ] , (A.3) 



1uj n R 

n=—oo 



where 

[An, A m ] = 5 n:rn , 

and 

j — — 27T??< . . 

u n = ypi + m 1 , p n = n = 0,±l,... (A.4) 

Using the above expansion together with the commutation relation of A and A\ it is easy 
to compute the propagator of the field, given by 

oo 1 

A F (x-x',t-t') = (<P{x,t)(f){x',t')) = V -iK(*-*')-M*-*')] . (A.5) 

' 2,Lu n R 

n=— oo 

The vacuum expectation value of the operator 2 (x, t) is then formally given by 

(0 2 (x,t)) = A F (0) (A.6) 

and, by translation invariance, is independent from x and t. However this expression 
is divergent and needs therefore to be regularized. Analogously to what has been done 
in the text for the ground state energy <?Q ac (.R), we need to subtract the corresponding 
expression in the infinite volume, so that the finite quantity, simply denoted by (f>o(R), 
satisfies the usual normalization condition 

hm ^ (R) = . 

H— >oo 
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Hence we define 



oo 



Isolating its zero mode, the series needs just one subtraction, i.e. 

oo oo 1 1 I oo 1 

71=1 



n=l y^T^j 1 n=l ( y^T^j 1 



(r = mR). In the above expression, the first series is now convergent whereas the second 
series, which is divergent, has to be combined with a divergence coming from the integral. 
Indeed we have 



X(r) = [ dn — 1 = lim (ln2A-ln — )- limlnA+ lim In A , (A. 9) 

J / 2 , ( r \ 2 A^oo I 2lX ) A^oo A^oo 



and the last term can be used to compose (|2.25)l . Collecting the above expressions, it 
is now easy to see that <Pq{R) coincides with the one obtained doing the calculation in 
the other channel, i.e. at a finite temperature. In fact, using the results of Ref. [21], this 
quantity can be expressed as 



a 27re rcosh0 - 1 



whose expansion in r is given by 



6KR) = _ + — flog — + 7g -log2 > ) f — 1 - — ) • (A.ll) 

Also this result could have been directly obtained computing only the finite part of 
the integral and using the prescription ([2.30j) . 



B Elliptic integrals and Jacobi's elliptic functions 

In this appendix we collect the definitions and basic properties of the elliptic integrals 
and functions used in the text. Exhaustive details can be found in [5Uj . 

The complete elliptic integrals of the first and second kind, respectively, are defined 

as 

tt/2 tt/2 

K(k 2 ) = f — da =, E(k 2 ) = [ daVl-k 2 sin 2 a. (B.l) 
J v 1 — k 2 sin 2 a J 
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The parameter k, called elliptic modulus, has to be bounded by k 2 < 1. It turns out 
that the elliptic integrals are nothing but specific hypergeometric functions, which can be 
easily expanded for small k: 

n 2 

k 2n + . . 



K(k 2 ) 
E(k 2 ) 



7T 



1 1 



7T 



1 1 '■> k 
2'2' ' 



■ « « 1 , k 
2'2' 



7T 

2 



l + ifc 2 + -A; 4 
4 64 



(2n-l)!! 



2 14 64 



2 n n! 

(2n-l)!! 
2™nJ 



2 fc 2n 



2n - 1 



Furthermore, for k 2 —> 1, they admit the following expansion in the so-called complemen- 
tary modulus k' = \/l — k 2 : 



K(k 2 ) 
E(k 2 ) 



log-+ (log- 



k> 



.12 



4 



+ 



(B-3) 



4 1\ k' 2 

|log F-2 T + 



Note that the complementary elliptic integral of the first kind is defined as 

K'(k 2 ) = K(k' 2 ) . 



(B-4) 



The function am(u, k 2 ), depending on the parameter k, and called Jacobi's elliptic 
amplitude, is defined through the first order differential equation 

2 

= 1 - k 2 sin 2 [am(u)] , (B.5) 



da,m(u) 



du 

and it is doubly quasi-periodic in the variable u: 

am (u + 2n~K + 2zmK') = rm + am(w) . 
The Jacobi's elliptic function sn(rt, /c 2 ), defined through the equation 

is related to the amplitude by snu = sin (amw), and it is doubly periodic: 

sn (u + 4nK + 2zmK') = sn(u) . 



(B.6) 



C Lame equation 

The second order differential equation 

^-E-N(N+l)V(u)}f(u) = 0, (C.l) 
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where E is a real quantity, N is a positive integer and V{u) denotes the Weierstrass 
function, is known under the name of iV-th Lame equation. The function V(u) is a 
doubly periodic solution of the first order equation (see [371] ) 



du 



4(V-e 1 ) (V-e 2 ) (V - e 3 ) , (C.2) 



whose characteristic roots e±, e2, uniquely determine the half-periods u> and u', defined 
by 

V (u + 2racj + 2mu') = V{u) . 

The stability equation ()3.7|) can be identified with eq. (jC.lj) for N = 1, u = x+iK.' and 
E = ^y- — k 2 cu 2 in virtue of the relation between V{u) and the Jacobi elliptic function 
sn(w, k) (see formulas 8.151 and 8.169 of [371]): 

fc 2 sn 2 (x, jfe) = V(x + iK') + . (C.3) 

Relation ()(I3[) holds if the characteristic roots of V(u) are expressed in terms of k 2 as 

2-k 2 2k 2 -1 1 + k 2 
ei = — g — , e 2 = — , e 3 = — , (C.4) 

and, as a consequence, the real and imaginary half periods of V(u) are given by the elliptic 
integrals of the first kind 

u = K(Jfe) , J = iK'(k) . (C.5) 

All the properties of Weierstrass functions that we will use in the following are specified 
to the case when this identification holds. 

In the case N = 1 the two linearly independent solutions of ()C1|) are given by (see 

EH) 

ha(u) = e=F«CW , (C.6) 

a(u) 

where a is an auxiliary parameter defined through V(a) = E, and o(u) and ((u) are other 
kinds of Weierstrass functions: 

*M = -V(u) , li^M = , (C.7) 

with the properties 

C(n + 2K) = C(«) + 2C(K) , 

a{u + 2K) = - e 2(u+K)f(K V(n) . (C.8) 
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As a consequence of eq. ()C.8|) one obtains the Floquet exponent of f± a (u), defined as 



/(u + 2K) = f(u)e iF ^, (C.9) 

in the form 

F(±a) = ±2i[KC(a)-aC(K)] . (CIO) 

The spectrum in the variable E of eq. (|C.1J) with N = 1 is divided in allowed/forbidden 
bands depending on whether F(a) is real or complex for the corresponding values of a. 
We have that E < e% and €2 < E < e\ correspond to allowed bands, while e% < E < 
and E > e\ are forbidden bands. Note that if we exploit the periodicity of V(a) and 
redefine a — > a' = a + 2nu + 2mu>' , this only shifts F to F' = F + 2m7r. 

The function admits a series representation [32] that will be very useful for our 
purposes in Sect. 



. 71 /"ku\ /E k 2 - 2\ 27r ^ h 2n . frmu\ 

C(u) = cot + 1 ) u H > — sm , (C.ll) 

SV ; 2K V2K/ VK 3 J K^l-h 2n VK/' V ; 

x ' n=l 

where h = e _7rK '/ K . The small-A; expansion of this expression gives 

£(u) = (cotu + -J H (u - 3cotu + 3ucot 2 u) + (C.12) 

V 3 / 12 

+ — (— 3u + (4u 2 — 5) cot u + u cot 2 u + Au 2 cot 3 it + sin 2uj + . . . 

(note that h ~ (|) 2 + O (k 4 )). A similar expression takes place for V(u), by noting that 
V{u) = - d -^. 
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